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Abstract 

We measure the inclusive branching fractions of charm mesons into three mesons with large sis 
content, namely the n, n' and 4>. Data were accumulated with the CLEO-c detector. For D° and 
D + rates, we use 281 pb _1 taken on the ^(3770) resonance, and for Df rates, we use 195 pb _1 
taken at 4170 MeV. We find that the production rates of these particles are larger in Df decays 
than in D° and D + decays. The cj> rate, in particular, is 15 times greater. These branching fractions 
can be used to measure B s yields either at the T(55) resonance or at hadron colliders. 
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I. INTRODUCTION 

Inclusive decay rates of charm mesons into mesons with large ss inherent quark content, 
the 77, rf and <ft mesons, are important for studies of both charm and b decays. Non-strange 
D mesons are generally expected not to decay into such objects, while the is likely to 
have both s and s quarks present after the primary c — > sW + transition, resulting in many 
more such particles. This is particularly useful in distinguishing between B and B s mesons 
as the B decays into D's with a large rate of ~90%, while the decays into are at the 
~10% level. We expect that the reverse is true for B s mesons. Knowledge of the charm 
yields into these mesons would allow alternative analyses of B s rates at the T(55) or at 
hadron colliders pj. 

In this analysis we use 281 pb -1 integrated luminosity of CLEO-c data produced in e + e~ 
collisions and recorded at the peak of the ip" resonance (3.770 GeV) to study the 77, rj' and 
4> yields in D° and D + decays. Production in decays is studied at 4170 MeV, where the 
cross-section for D*^Df is ~1 nb 0. 

The CLEO-c detector is equipped to measure the momenta and directions of charged 
particles, identify charged hadrons, detect photons, and determine with good precision their 
directions and energies. It has been described in more detail previously [3j, \% |5( . 



II. SELECTION OF D°, D~ AND Dj TAGGING MODES 

Fully reconstructed charged or neutral D meson candidates are selected from the data at 
3.770 GeV, where pairs of D°D or D+D- mesons are produced. The decay modes used 
are listed in Table |1J In general, in this paper, D refers to either a D° or D + meson and 
its anti-particle, and D s refers to meson and its anti-particle. (Also, mention of a flavor 
specific state also implies use of the charge-conjugate state.) 



Mode 


Signal 


Background 




77387 ± 281 


1868 




24850 ± 214 


12825 


K S 7T~ 


11162 ± 136 


514 




18176 ± 255 


8976 


K S 7T-TT° 


20244 ± 170 


5223 


Sum D~ 


151819 ±487 


29406 


K+TT~ 


49418 ± 246 


630 


K+1T-TT 


101960 ± 476 


18307 




76178 ± 306 


6421 


Sum 


227556 ± 617 


25357 



TABLE I: Tagging modes and numbers of signal and background events determined from the fits to 
the D~ and D beam constrained mass distributions, after making the mode-dependant tobc cuts. 
The error on the summed signal yield is obtained by adding the errors on the individual yields in 
quadrature. 
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At 4.170 GeV we produce DfD~ pairs, with one of the D s being, most of the time, the 
daughter of a D* decay. Fully reconstructed Df candidates are selected in the decay modes 
listed in Table |H] D mesons at this energy are a source of background, they are mostly 
produced in D*D final states, with a cross-section of ~5 nb, and D*D + DD final states, 
with a cross-section of ~2 nb. DD is a relatively small, ~0.2 nb. There also appears to be 
DD production with extra pions. 



Mode 


Signal 


Background 


K+K--K- 


8446 ± 160 


6793 


K S K- (K s -» 7T+7T-) 


1852 ± 62 


1022 




1101 ±80 


2803 




786 ± 37 


242 


4>p- 


1140 ± 59 


1515 


K*+(890)K*°(890) 


1197 ±81 


2599 


Sum 


14522 ± 218 


15328 



TABLE II: Tagging modes and number of signal and background events determined from the fits 
shown in Fig. |1J after making the mode dependant invariant mass cuts. The error on the summed 
signal yield is obtained by adding the errors on the individual yields in quadrature. 

We fully reconstruct one of the D mesons at 3.770 GeV or one of the D s mesons at 4.170 
GeV to form a specific tag, and then look for cases where the particle produced in association 
with our tag has a decay of either rj — > 77, 7/ — > 7i + 7i^i], i] — > 77, or — > K + K~ . 

All acceptable track candidates must have a helical trajectory that approaches the event 
origin within a distance of 5 mm in the azimuthal projection and 5 cm in the polar view, 
where the azimuthal projection is in the bend view of the solenoidal magnet. Each track 
must possess at least 50% of the hits expected to be on a track, and it must be within the 
fiducial volume of the drift chambers, | cos#| < 0.93, where 6 is the polar angle with respect 
to the beam direction. 

We reconstruct 7r°'s by first selecting photon candidates from energy deposits in the 
crystals that are not matched to charged tracks and that have deposition patterns consistent 
with that expected for electromagnetic showers. Pairs of photon candidates are kinematically 
fit to the known 7r° mass 0]. We require the pull, the difference between the reconstructed 
and known tt° mass normalized by its uncertainty, to be less than three for acceptable 7r° 
candidates. 

Ks candidates are formed from a pair of charged pions that are constrained to come from 
a single vertex. We also require that the invariant mass of the two pions be within 4.5 times 
the width of the Ks mass peak, which has an r.m.s. width of 4 MeV. 

We use both charged particle ionization loss in the drift chamber (dE/dx) and Ring- 
Imaging Cherenkov (RICH) information to identify kaons and pions used to fully reconstruct 
D and D s mesons. The RICH is used for momenta above 0.7 GeV/c. The angle of detected 
Cherenkov photons that were radiated by a particular charged track are translated into an 
overall likelihood denoted by Li for each particle hypothesis. To differentiate between pion 
and kaon candidates, we require the difference —2 log(£ 7r ) + 2 log(£^) to be less than zero. 
To utilize the dE/dx information, we calculate a n as the difference between the expected 
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ionization loss for a pion and the measured loss divided by the measurement error. Similarly, 
Ok is defined using the expected ionization for a kaon. 

We use both the RICH and dE/dx information for D and D s meson tag candidate tracks 
in the following manner: (a) If neither the RICH nor dE/dx information is available, then the 
track is accepted as both a pion and a kaon candidate, (b) If dE/dx is available and RICH is 
not then we insist that pion candidates have PIDje = cr^— c K < 0, and kaon candidates have 
PIDdE > 0. (c) If RICH information is available and dE/dx is not available, then we require 
that PID R ich = — 21og(£ 7r ) + 21og(£x) < for pions and PID RIC h > for kaons. (d) If 
both dE/dx and RICH information are available, we require that [PID^e + PI D rich) < 
for pions and [PlD dE + PID RIC h) > for kaons. 

A. Reconstruction of D and D~ Tagging Modes 

Tagging modes for D° and D + decays are reconstructed as described previously 0]. 
Briefly, at the ?/>(3770) D meson final states are reconstructed by first evaluating the differ- 
ence, AE, between the energy of the decay products and the beam energy. We fit the AE 
spectrum with a double Gaussian to represent the signal, and a polynomial representing the 
background. We require the absolute value of this difference to contain 98.8% of the signal 
events, i. e. to be within ~2.5 times the r.m.s width of the peak value. For final states 
consisting entirely of tracks, the AE resolution is ~ 7 MeV. A tt° in the final state degrades 
this resolution by roughly a factor of two. Candidates with AE consistent with zero are 
selected, and then the D beam-constrained mass is evaluated, 



where i runs over all the final state particles. 

The wbc distributions for all D and D~ tagging modes are shown in Figs. ^ and El 
Table lists the numbers of signal and background events within the signal region defined as 
containing 98.8% of the signal events with tubc below the peak and 95.5% of the signal events 
above the peak; the interval varies from mode to mode. The numbers of tagged events are 
determined from fits of the wbc distributions to a signal function plus a background shape. 
The signal is described by a Crystal Ball Line shape j^, • For the background, we fit with 
a shape function analogous to one first used by the ARGUS collaboration [8], which has 
approximately the correct threshold behavior at large wbc, except for the D° — > K~tt + and 
D° — > K~7t + 7t° modes where we use a fourth order polynomial. For each tagging mode, the 
background function is first fit to a tti-qq distribution that lies within an interval from 5 to 
7.5 r.m.s. widths away from the peak of the AE distribution. We fix the shape parameters 
from these fits and then use these functions for background distributions in the signal fits, 
allowing the normalization to float. 

We find 151819 ± 487 ± 759 D~ and 227556 ± 617 ± 1138 D° signal events that we use 
for further analysis. The systematic uncertainty on this number is estimated to be ±0.5% 
by varying the fitting functions. 
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FIG. 1: Beam-constrained mass distributions for fully reconstructed D° decay candidates in the 
final states: (a) K+vr", (b) K+ 7r 7r , and (c) K + ir it ir + . The distributions are fit to a Crystal 
Ball Line shape for the signal. For the background, we either use a fourth order polynomial (in 
(a) and (b)) or an ARGUS shape (in (c)). Both background shapes are obtained from the AE 
sidebands. The regions between the arrows are selected for further analysis. 

B. Reconstruction of Dj Tagging Modes 

At 4170 MeV the presence of the 7 from D*~ — > causes us to adopt a different 

procedure. If we ignore the photon and reconstruct the tti-qq distribution, we obtain the 
distribution from Monte Carlo shown in Fig. 01 The narrow peak occurs when the recon- 
structed D s does not come from the D* decay. Thus, the method of applying narrow cuts 
on mbc and AE, used so successfully on the ^(3770), no longer works. 

Instead, we insist that the D~ candidate has momenta which satisfies the requirement 
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FIG. 2: Beam-constrained mass distributions for fully reconstructed D~ decay candidates in the 
final states: (a) K + it~-k~ , (b) K + tt~tt~it®, (c) K s it~ , (d) K s it~tt + tt~ , and (e) K s tt~tt°. The 
distributions are fit to a Crystal Ball Line shape for the signal and an ARGUS shape obtained 
from the AE sidebands for the background. The regions between the arrows is selected for further 
analysis. j 




2.00 2.05 

m BC (GeV) 



2.10 



FIG. 3: The beam constrained mass tobc from Monte Carlo simulation of e + e~ — > D+D*~, 
D*~ — > 7-Dj, Djr — > (f)^ at 4170 MeV. The narrow peak is from the D£ and the wider one from 
the D~ . (The distributions are not centered at the or masses, because the reconstructed 
particles are assumed to have the energy of the beam.) 



2.015 < mec < 2.067 GeV. This requirement allows for the fact that the D~ could have 
been produced directly or as a result of a D* s ~ decay to either jDj, or tt°D~ decay with 
a small ~ 5.8% branching fraction jfj. We then reconstruct the invariant mass of the D s 
candidates. The invariant mass distributions of the tagging modes we considered in this 
analysis are shown in Fig. they are listed in Table |H1 where we also give the number of 
signal and background events. Here we use only the secondary decays i^*°(890) — > K + tt~ , 
K* + (890) — > K s tt + , p~ — > 7T~7r°, rj — > 77 and 7/ — > tt + 7t~7]. (More specifically, when 
appropriate, we require the Kir invariant mass to be within ±100 MeV of the K* mass, the 
7T~7r° mass to be within ±100 MeV of the p~ mass, the 77 invariant mass minus the 7] mass 
divided by its error to be less than 3, and the invariant mass of the 777r + 7r~ minus the 77 
mass, for 77 combinations consistent with the 7] hypothesis, to be within ±10 MeV of the 
known if — rj mass difference.) 

The D~ signal regions are defined as containing 98.8% of the signal events within an 
invariant mass window symmetric about the D~ mass peak. The intervals vary from mode 
to mode. To find the numbers of signal tag events, the invariant mass distributions are 
fit to a sum of two Gaussian signal functions with the means constrained to be the same 
and the r.m.s. widths allowed to float. The background function is a second or third order 
Chebychev polynomial. 

We have 14522 ± 218 ± 145 D~ signal events that we use for further analysis. The 
systematic uncertainty is estimated to be ±1.0% by varying the signal and background 
fitting functions. 
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FIG. 4: Invariant mass distributions for fully reconstructed D s decay candidates in the final states: 
(a) K+K-ir-, (b) K a K~, (c) 7771-", (d) t/tT, (e) <f>p~, and (f) K*°K*+. The distributions are fit to 
double Gaussian signal shapes and Chebychev polynomial backgrounds. The regions between the 
arrows is selected for further analysis. 
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III. 77, rf AND SELECTION 



For the 77 we use the 77 final state, which has a measured branching fraction of (39.43 ± 
0.26)%. To detect 7/ we use the 7r + 7r~77 final state, which has a branching fraction of (44.3 ± 
1.5)%, with the 77 subsequently decaying into 77. For the we use the K + K~ final state 
with a rate of (49.1 ± 0.6)% @. 

The track selection and particle identification requirements for the signal side are identical 
to those for the tag side, except for momenta less than 0.2 GeV/c, where we loosen the dE/dx 
consistency requirement to 4 ax- This is the case for both D and D s meson tags. 

We accept photons only in the best-resolution region of the detector, | cos#| < 0.8, where 
9 is the angle of the photon with respect to the beam direction. Photon candidates must not 
be matched to charged tracks, must have a reconstructed energy greater than 30 MeV and 
have a spatial distribution in the crystals consistent with that of an electromagnetic shower. 

Candidates for 77' mesons are selected by combining 77 candidates within 3 r.m.s. widths 
of the 77 mass, with a tt + and a 7r~. The mass difference between f]ir + TT~ and 77 is then 
examined. Both pions forming 7/ and kaons forming candidates are required to pass the 
track selection and particle identification requirements. The signal yield is then extracted 
from fits to the 77, 7/ — 77, and mass distributions. 

IV. RECONSTRUCTION EFFICIENCIES 

The reconstruction efficiencies for 77, 77' and in our tag samples of D and D s events are 
shown in Fig. |3| They are determined from a Monte Carlo simulation of the detector [ill ]. 
There is no observable difference in the efficiencies for D° and D + decays. In the case of the 
D s , the reconstruction efficiency, especially for the 0, is lower in the highest momentum bin, 
due to the different angular distributions of these particles, caused by the different production 
angles with respect to the beam of charmed mesons from DD production compared with 
D*D S production. In other words, the kaons on average are produced closer to the beam 
axis in D*D S events than in DD events. 

The 77 efficiency falls slowly below 300 MeV/ c, and then levels out. Since our aim here is to 
measure the inclusive branching fractions, we break the 77 sample into two parts, one below 
300 MeV/c and the other above. For the 77', the efficiency is constant with momentum, 
so we do not separate the data into momentum intervals. The efficiency, on the other 
hand, changes drastically with momentum and therefore we use several momentum regions. 
The increase in the efficiency is easily explained by the fact that as the becomes more 
energetic it is less likely to produce a kaon of p < 0.2 GeV/c, which would cause the event 
to be rejected. 

The simulated efficiency could be inaccurate if the polarization were not correct; this 
could occur because of a poor choice of the mixture of final states. The data and Monte Carlo, 
however, show the same polarization. (The observed polarization is almost independent of 
momentum in both data and Monte Carlo.) 
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FIG. 5: Reconstruction efficiencies for: r/ — > 77, 7/ — > r/7r + 7r~, and (/> — > i"T + i<C~. The filled circles 
indicate D°D° events, the open circles D + D~ events, and the crosses D*D S events. The efficiencies 
do not include branching ratios. 

V. SIGNAL YIELDS AND BRANCHING FRACTIONS 



The signal yields in this analysis are evaluated by taking the difference between the 77, rf 
and yields opposite selected tags and the yields in samples that estimate the background 
under the tag peaks. Our procedure is somewhat different for D and D s decays. In the D 
case, we evaluate the background yields using events in the low and high sidebands of the 
AE distribution from 5a to 7.5a away from the peak. These sidebands are normalized to 
have the same number of events as the backgrounds under the AE peaks. In the D s case we 
select sidebands in the same interval relative to the peak as for the D but in invariant mass 
rather than in AE. 



A. Inclusive 77 Yields 

In Fig. El we show the two-photon invariant mass in our two momentum intervals for both 
signal and sideband regions for D tags; Fig. [7| shows the corresponding distributions for D s 
tags. 

The 77 signal and background yields are determined by fits to a Crystal Ball function, to 
account for the peak and the low mass tail, and a background polynomial. For the signal 
D° region, the four fit parameters describing the Crystal Ball line shape, three for shape 
and one for the mean are allowed to float. These parameters are then fixed and used for the 
D° sideband regions and the D + and signal and sideband regions, since some of these 
have relatively small yields. Table IIHI lists the yields, the efficiencies, and the branching 
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FIG. 6: Invariant mass of the 77 — > 77 candidates from D° decay: (a) signal region events with the 
momentum of 77, \p v \, less than 0.3 GeV/c, (b) signal region events with 0.3 < \p v \ < 1.0 GeV/c, 
(c) sideband events with \p v \ < 0.3 GeV/c, (d) sideband events with 0.3 < \p v \ < 1.0 GeV/c. 
Candidates from D + decay are shown in (e)-(h), with corresponding descriptions. 



fractions for the two momentum intervals. (Yields in the highest momentum bin include a 
small number of events that slightly exceed 1 GeV/c.) 

The systematic uncertainties arise from several sources. For the 77 we estimate a detection 
efficiency error of ±2% per photon 1 for a total of ±4%. We also include an error due to fitting 



1 This is determined from a study of pions in ip(2S) — > J/ipir°ir° transitions |12(. 
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FIG. 7: Invariant mass of the r\ — > 77 candidates from decay: (a) signal region events with the 
momentum of 77, \p v \, less than 0.3 GeV/c, (b) signal region events with 0.3 < \p v \ < 1.0 GeV/c, 
(c) sideband events with \p v \ < 0.3 GeV/c, (d) sideband events with 0.3 < \p v \ < 1.0 GeV/c. 



the Monte Carlo samples. In addition, there is an uncertainty caused by using the efficiency in 
only two momentum intervals, due to possible variations in these intervals, amounting to 3%. 
This error is estimated by considering the effects of different parent momentum distributions. 
The uncertainties on the tag yields are derived by varying the fitting functions. We also have 
significant contributions from uncertainties on the signal and background fitting function, 
determined by varying the functions. All the systematic error contributions are listed in 
Table IIVI For the 7/ and 4> modes we also list the estimated uncertainties for finding the 
charged tracks and identifying their species. These differ somewhat between the rj' and 
because of the different track momenta involved. For the mode we include another 
additional error source due the lack of efficiency in the first momentum bin that we will 
discuss in more detail subsequently. 
For the 77 rates we find 

B(D° -> r]X) = (9.5 ± 0.4 ± 0.8)% (2) 
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Tag 1 


In 1 
\Pv\ 


M S1 S 


/V bk S 
Jv »7 


TV e % 


BHD — >■ nX) 




(GeV/c) 








(%) 


n° 

i j 


on n ^ 

U.U - U.O 


1454 ± 133 


176 ± 33 


1278 ± 137 57 


O K -L Q 

Zi.o in u.o 




0.3-1.0 


3427 ± 137 


242 ± 36 


3185 ± 141 50 


7.0 ± 0.3 


Sum 




4880 ± 191 


418 ±49 


4463 ± 197 - 


9.5 ±0.4 


D+ 


0.0-0.3 


608 ± 65 


153 ± 35 


455 ±74 58 


1.3 ±0.2 




0.3-1.0 


1811 ± 115 


294 ± 39 


1517 ±121 51 


5.0 ±0.4 


Sum 




2419 ± 132 448 ± 332 1972 ± 142 - 


6.3 ±0.5 


Df 


0.0-0.3 


230 ± 38 


152 ± 36 


78 ± 53 55 


2.5 ± 1.7 




0.3-1.0 


963 ± 56 


367 ± 48 


596±74 50 


21.0 ±2.6 


Sum 




1193 ± 68 


519 ±60 


674 ±91 - 


23.5 ±3.1 



TABLE III: r] signal yields (A^ lg ), background yields (iV^ kg ) and background-subtracted yields (N v ) 
as a function of momentum. Also listed are the rj reconstruction efficiencies (e l ) in percent, and 
the partial branching fractions versus momentum. (Yields in the highest momentum bin include a 
small number of events that slightly exceed 1 GeV/c.) 



Systematic uncertainties 


V (%) 


rf(%) 


<K%) 


Photon reconstruction 


4.0 


4.0 




Charged track finding 




1.4 


5.0 


Particle identification 




2.0 


2.0 


Monte Carlo fitting 


2.0 


1.0 


1.0 


Average efficiency 


3.0 


3.2 




Number of tags (D° & D+) 


0.5 


0.5 


0.5 


Number of tags (Df) 


1.0 


1.0 


1.0 


Signal & Background Fitting 


6.5 


5.3 


2.1 


Estimate of 1st p bin (D° k D + ) 






2.0 


Estimate of 1st p bin (Df) 






3.1 


Total (D° & D + ) 


8.5 


8.6 


6.3 


Total (D+) 


8.5 


8.6 


6.8 



TABLE IV: Systematic uncertainties (±) on the inclusive rj, rj and (ft branching ratios. 

B(D + ^r]X) = (6.3 ±0.5 ±0.5)% 
B(Df^rjX) = (23.5 ±3.1 ±2.0)%. 

Note that these rates naturally include cascade decays from rj' — > r/X. 
B. Inclusive r/ Yields 

We first reconstruct 77 invariant mass as shown above. Then we use 77 mass combinations 
within ±3<T of the rj mass as 77 candidates, where a is the r.m.s. width of the mass peak. In 
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Figs.|Hland|niwe show the i]7t + tt~ - r\ mass difference for D and D s tags from both signal and 
sideband regions. To determine the yields in this case we fit to a Gaussian signal function 
and a background polynomial. The signal shapes and means are allowed to float for the 
signal distributions and fixed to the values obtained there in the corresponding sideband 



regions. 
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FIG. 8: Difference in the invariant mass of rf — » r]ir + Tr~ and rj (rj — > 77) candidates from: (a) 
D° signal region events (b) D + signal region events (c) D° sideband events, and (d) D + sideband 
events. The fits are described in the text. 



The signal, background, and background-subtracted yields, the detection efficiency and 
the branching fraction are given in Table The systematic error sources are listed in 
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FIG. 9: r]ir + Tr~ - rj mass difference from DfD~ signal events (a) and sideband events (b). The 
fits are described in the text. 

Table HVI For the charged tracks we estimate a systematic error of ±0.7% for track finding 2 
and ±1% for particle identification. 



Tag 


N m , g 

V 




N v , 


e* 


B(D (s) ^r,'X)(%) 


D° 


313 ± 20 


14 ±5 


299 ± 21 


30 


2.48 ±0.17 ±0.21 


D+ 


83 ± 12 


lit 


82 ± 13 


30 


1.04 ±0.16 ±0.09 


Df 


91 ± 12 


23 ±8 


68 ± 15 


31 


8.7 ± 1.9 ±0.8 



TABLE V: rf signal yields (N^, g ), background yields (N^, g ) and background-subtracted yields 
{Nrf), the rj' reconstruction efficiencies (e*), and the branching fractions. 



C. Inclusive Yields 

In Figs. [TUllTol we show the K + K~ invariant mass for the signal region in five different 
momentum intervals from both signal and sideband regions for D°, D + and tags respec- 
tively. The signals are fit with a sum of two Gaussian shapes and the background is fit to 
a polynomial. The signal shapes are fixed to the values obtained by Monte Carlo simula- 
tion, while the mean is allowed to float. The signal, background, and background-subtracted 



2 This is determined from a study of pions in ijj{2S) — » J/ipir + ir transitions [12| but increased to account 
for the contribution from low momentum tracks. 
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1 ns + n 1 9 
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0.0 - 


0.2 


1.4 it 1.6 


1.2 + 1.4 


0.1 + 2.1 


3.7 


0.1 + 0.8 


0.2 - 


0.4 


49.0 ±7.4 


8.3 ±3.5 


40.7 + 8.2 


14.6 


3.9 + 0.8 


0.4 - 


0.6 


144.4 ± 12.9 


28.4 ±6.2 


116.1 ± 14.3 


32.1 


5.1 + 0.6 


0.6 - 


0.8 


155.6 ± 13.1 


20.3 ±5.1 


135.3 + 14.1 


44.7 


4.2 + 0.5 


>0.8 




82.1 ±9.1 


6.3 ±3.0 


75.8 ±9.6 


44.4 


2.4 + 0.3 


Total 


432.5 ±21.9 


64.5 ±9.4 


368.0 ± 23.8 




15.7 + 1.4 



TABLE VI: Measured 4> signal yields (N^ g ), background yields (N^ g ) and background-subtracted 
yields (JVa) versus momentum from D and D s decays. Also listed are the (ft reconstruction efficien- 
cies (e 4 ), and the partial branching fractions vs momentum. (Note that measurements in the lowest 
momentum interval will be replaced by a model dependent estimate.) 

yields, the detection efficiency and the branching fraction in each momentum interval are 
given m Table EH For D° and D + there are no events above 0.9 GeV/c, while for the -D+ 
there small number of events above 1.0 GeV/c. 

Although the measured yields in the lowest momentum bin, 0< p <0.2 GeV/c, are quite 
small, so are the efficiencies. To take into account possible incorrect efficiency estimates in 
this difficult kinematic region we make an independent estimate of the rate by using the 
Monte Carlo predicted fraction of the <fi yield. We then take a conservative 100% error 
on these estimates. The fractions are 2.0%, 2.0% and 3.1%, for D°, D + and D+, respec- 
tively. These correspond to partial branching fractions of (0.02±0.02)%, (0.02±0.02)%, and 
(0.5±0.5)%, respectively. Using these more reliable estimates for the rates in the first bin, 
we show the efficiency corrected momentum distributions in Fig. The inclusive branching 
ratios are 

B(D°^(f)X) = (1.05 ±0.08 ±0.07)% (3) 
B(D + ^(f)X) = (1.03 ±0.10 ±0.07)% 
B(Df^(f)X) = (16.1 ± 1.2 ± 1.1)% . 
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The systematic uncertainties are listed in Table II VI Due to the difficulty in tracking slow 
charged kaons, we use an increased systemic error of ±2.5% per track in evaluating the 
efficiency error from track finding. The error due to particle identification remains at ±1%. 

VI. CONCLUSIONS 

Our results of the inclusive rj, rj' and <p production rates from D°, D + and £)+ decays are 
summarized in Table IVHI Of the 9 measured rates in this paper, 8 are first measurements 
and the other one, D° <f)X, improves the accuracy from 50% to 10%. We are consistent 
with previous upper limits |6[ in the three cases where they exist. 



Mode 


D°{%) 




D+{%) 




£>+(%) 




Our result 


PDG 


Our result 


PDG 


Our result PDG 


rj X 


9.5±0.4±0.8 


<13 


6.3±0.5±0.5 


<13 


23.5±3.1±2.0 


7]' X 


2.48±0.17±0.21 




1.04±0.16±0.09 




8.7±1.9±0.8 


(f> X 


1.05±0.08±0.07 


1.7±0.8 


1.03±0.10±0.07 


<1.8 


16.1±1.2±1.1 



TABLE VII: Summary of inclusive branching ratio results. 

These particles all have significant components of ss. Our results show that rj' and <p are 
relatively rare in D° and D + decay while the rj which has a lower mass and a significant 
light quark component, is produced at a significantly higher rate. The rj, rj' and <fi are all 
produced at higher rates in D s decays than the corresponding rates from D decays. The 
ratio of rates is given in Table IVTiTl The 4> yield is 15 times higher in decays than in D 
decays. 



Ratio 


V 


rf 


<t> 


Df/D° 


2.47 ±0.34 ±0.18 


3.51 ±0.80 ±0.27 


15.3 ± 1.6 ±0.8 


D+/D+ 


3.73 ±0.57 ±0.27 


8.37 ±2.23 ±0.64 


15.6 ± 1.9 ±0.8 



TABLE VIII: Ratios of Dj to D yields. Common systematic errors have been eliminated. 

The large asymmetry in the yields of these particles between D s and the lighter D mesons 
will permit further studies of B s decays at the T(5S), and will be most useful for separating 
B s from B decays at hadron colliders. 
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FIG. 10: Invariant mass of <j> — > K + K candidates from D°D° signal events in five different 
momentum intervals: (a) < < 0.2 GeV/c, (b) 0.2 < < 0.4 GeV/c, (c) 0.4 < < 
0.6 GeV/c, (d) 0.6 < \p^\ < 0.8 GeV/c, (e) 0.8 < |p^| < 0.9 GeV/c. The fits are described in the 
text. 
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FIG. 11: Invariant mass of (ft — > K + K~ candidates from D°D° background events in five different 
momentum intervals: (a) < < 0.2 GeV/c, (b) 0.2 < < 0.4 GeV/c, (c) 0.4 < \p^\ < 
0.6 GeV/c, (d) 0.6 < {p^ < 0.8 GeV/c, (e) 0.8 < < 0.9 GeV/c. The fits are described in the 



text. 
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FIG. 12: Invariant mass of (f) — ► K + K candidates from D + D signal events in five different 
momentum intervals: (a) < \p^\ < 0.2 GeV/c, (b) 0.2 < \p$\ < 0.4 GeV/c, (c) 0.4 < \p^\ < 
0.6 GeV/c, (d) 0.6 < \p^\ < 0.8 GeV/c, (e) 0.8 < |p^| < 0.9 GeV/c. The fits are described in the 
text. 
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FIG. 13: Invariant mass of <f> — > K + K candidates from D + D background events in five different 
momentum intervals: (a) < \p^\ < 0.2 GeV/c, (b) 0.2 < \p^,\ < 0.4 GeV/c, (c) 0.4 < < 
0.6 GeV/c, (d) 0.6 < \p^\ < 0.8 GeV/c, (e) 0.8 < |p^| < 0.9 GeV/c. The fits are described in the 
text. 
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FIG. 14: Invariant mass of 4> — > K + K~ candidates from DfDj signal events in five different 
momentum intervals: (a) < \p^\ < 0.2 GeV/c, (b) 0.2 < \p^\ < 0.4 GeV/c, (c) 0.4 < \p^\ < 
0.6 GeV/c, (d) 0.6 < |p^| < 0.8 GeV/c, (e) \p$\ > 0.8 GeV/c. The fits are described in the text. 
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FIG. 15: Invariant mass of <f> — > K + K candidates from DfD s background events in five different 
momentum intervals: (a) < < 0.2 GeV/c, (b) 0.2 < < 0.4 GeV/c, (c) 0.4 < |p^[ < 
0.6 GeV/c, (d) 0.6 < \p^\ < 0.8 GeV/c, (e) 0.8 < |p | GeV/c. The fits are described in the text. 
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FIG. 16: The branching fraction in % in each 200 MeV interval of momentum for <j) mesons from 
(a) D° decays, (b) D + decays and (c) decays. 
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